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INTERACTION EET'JEEN THE SPANS OF SEMIMONOCOQUE WII'iGS WITH CUTOUTS 
By N. J. Eoff, Hari’y Kase, and Harold Llebowitz 

STOMAEY 


The stresses in the two spars of a model of a wing having three 
rectangular cutouts were calculated hy the PIB/iL method^ a modification 
of Southwell’s method of systematic relaxations. The model was built 
and tested in the Polytechnic Institute of Brooklyn Aeronautical 
Laboratories and the deflections and strains measured were compared 
with calculated values. The agreement was found to be satisfactory. 


INTRODUCTION 


Many modern airplane wings are built with two spars rather than ac- 
cording to the true monocoque principle. Spars are provided because the 
full utilization for load- -carrying purposes of the skin and its rein- 
forcements is not possible when there are many cutouts for retracting 
landing gear, wing tanks, armament, and so forth near the root. The 
calculation of the interaction between the spars is a difficult task if 
the rigorous methods of the theory of elasticity are used. To disregard 
it, however, is wasteful since the ribs and the skin provide a compar- 
atively strong elastic connection between the spars which relieves the 
more highly loaded spar at the expense of the less highly loaded one. 

The interaction between the spars of fabric-covered wooden >n.ngs was 
investigated by many authors. The earliest attempt appears to be con- 
tained in a paper by L. Ballenstedt preference 1) published in 1918^ 
while in the late twenties Th, von Karman’s group at Aachen developed 
and verified a rigorous and comparatively simple method of calculation 
(references 2 and 3). In England, D. Williams and E. Eoxbee Cox worked 
out an intex-esting solution in 1933 (reference 4) . All these papers 
dealt with wing structures without a load— carrying skin since at the 
time of their publication stressed- skin wings were seldom if ever used. 
However, in 1935 Paul Kuhn showed that the Friedrichs-Kannan equations 
can also be applied to the interaction problem of the more modern types 
of wing (reference 6) . More recently, a theoretical solution was given 
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"by V7. J. Goodey, who used strain-energy methods in his analysis (ref- 
erence 6) . 

The writers helieve that there is a need for a procedure which can 
take into account variahle cross sections, stressed skin, and cut- 
outs, and at the same time is simple enough to he used hy little— trained 
personnel. It appears that the PIBAL method developed in reference. ( 
for the calculation of the stresses in reinforced flat and curved panels, 
rings and frames, and reinforced monocoque cylinders is well suited for 
the solution of the interaction prohlem with an accuracy sufficient for 
engineering pui’poses. It is oa.sed on the Southwell relaxation mnthod 
(reference 8) . The structure is assuiaed to he composed of several heam 
elements having hending, shearing, and torsional rigidity determined 
from the geometry of the wing structure and tlie mechanical properties of 
its materials of construction. The forces and moments corresponding to 
prescribed displacements of the end points of these elements are deter- 
mined, and the conditions of equilihrium at these points are expressed 
in terms of the displacements. The result is a system of linear equa- 
tions which in the present report is solved hy matrix methods. In the 
Southwell method a solution is obtained hy a systematic procedure of 
step— hy— step approximations. The solution yields the displacements at 
the end points of the elements from which the stresses in any part of 
the structure can he easily calculated. 

The res-olts of the calculations were checked hy experiments carried 
out with a model of a two— -spar wing under various conditions of loading 
and end fixation. The agreement was found to he satisfactory. 

The authors are much indebted to Dr. Bruno A. Boley and Mr. 

Bertram Klein for their advice and help during the construction and 
testing of the specimen and in the calculations, to Edo Aircraft 
Comoration, and to Mr. R. Eies of Edo Aircraft Corporation for con- 
tributing the foimed rib flanges. The investigation was sponsored hy 
and conducted with the financial assistance of the National Advisory 
Committee for Aeronautics. 


SYlfflOIS 


A 


A 


c 


^i 



cross-sectional area of upper or lower flange 
cross-sectional area of compression flange 

area included hy the thin wall of a closed section 
cross-sectional area of uprights, of tension flange 
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Ax. A* 


area of cross section of an interspar element 


A£^ £A^ 33^ 
33, 35, etc, 


eutscripteused with influence coefficients, in vhich the 
first letter or nmaber refers to the point at which 
the force or moment exerted upon the "bar is acting 
eind the second to the point at which the displacement 
or rotation occurs 


h 

d 

E 


wall thickness 

deflection at end of beam element; distance between spar 
uprights (web stiffeners) 

Young’s modulus 

bending stress in cvirved sheet of front spar 


fB 

f^ 


maximum normal stress in element 3— 5 at its midpoint 


normal stress in compression flange due to f,^^ 


Cl 


normal stress in compression flange of front spar 
element due to bending resulting from T and Pi 

normal stress in compression fleinge due to bending 
resulting from T 

shear stress in curved sheet of front spar element 
normal stress in tension flange due to f^^^ 


ti 


normal stress in tension flange of front spar element 
due to bending resulting from T and Pi 

normal stress in tension flange due to bending resulting 
from T 


f,^ compressive stress in the uprights 

f^ tensile stress in the web parallel to the direction of 

wrinkles 

G shear modulus 

h distance between centroids of spar flanges 


4 


MCA TE^ rio. 1324 


2hi 


I 




I 


F 


I 


X 


tota.l length of uprights in a spar element 
moment of inertia 

moment of inertia of cui'ved sheet and flanges of, front 
spar 

moment of inertia of interspar element at front spar veh 

moment of inertia of intersrar element hetween spar wets 
at distance x forward of rear spar 



K 

L 


moment of inertia of interspar element between web and 
shear center of front spar, at distance x’ aft of 
shear center 

constant used in equations ( 9 ) 
length of beam element 


m 


rotation at end of beam element 


influence coefficients 


M 

Ms 

nn 


end moment acting on "beam element; bending noment at any 
section of front spar element due to T and Pi 

bending moment at the midpoint of element 3—5 


moment acting on element 3—5 ^^oint 3 
influence coefficient 


ni ^ n^^ 

n^ , n^., n^ rotations at ends of beam elements 

4 ' 5 -^ 6 


^2, Ns, 

II 4 , Ng, Ng end moments 


P 


total vertical load at end of front spar beam element 


P' 


perimeter of the cross section of an interspar element 
between the web 'and shear center of front spar 


Pi 


vertical load taken by curved sheet of front spar beam 
element 
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Ps 

vertical load on e.lement 3—5 Q't joint 3 


perimeter of the cross section of an inter spar element 
between spar webs, at distance x foDrwa.rd of rear 
spar 

Q 

static moment of area used in shear^-stress formula 
fg = PQ/It 

E ,H ,S , j ^ 
E.H.S.jj 

right-hand sides of equations of equilibrium, in which 
the subscripts x-efer to the condition of loading and 
end fixation 

s 

perimetric coordinate 

t 

web thickness 

t, ti, tg, 

CO 


t^, tg, tg rotations at ends of 1)68111 elements 
tt, ty influence coefficients 

'^St “^3} 



end moments 

0- 

\ total strain energy in a beam element 


^curved sheet strain energy in curved sheet of front spar elenient 


^flanges 

strain energy in flanges of front spar element 

^uprights 

strain energy in uprights of front spar element 

^web 

strain energy in flat web of front spar element 

V 

vertical displacement 

vv, vm 

influence coefficients 

V 

shear force; volume 

X, X* 

horizontal coordinates 
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vertical coorilnate 

vertical diaplace-uents 
influence coefficients 

Ye, 

sheai* forces 

i' 5 ' 6 

a angle of folds of spar webs 

0 ^ angular coci'di nates 

(0 rotation 

MCIHOD OF ANALYSIS 


y 

7} yi.» J2> 

y 4 > y 5 > Ys 
yy> 

Y, Yi, Ya, 
Y . Y . Y 


Fiaure 1 shows the wing structure that was analyzed ^d tested, exA 

figxires 2 and 3 are photographs of the wing on asSmLy 

noses of the ana-lysis the actual structure was replacea 
of bars shown in figui’e 4. The boxes consisting of two 
as well as the wing covering between them in tne a^ual 
replaced in this figure by bars 1-2, 3-^, 5 6. 

looted at the chcrdwise center lines of the bo-xes and 
shear center of the L-shape front spar to the Plane of f 

rear spar The spar elements extend spanwise from the cente. of on 
JS; centa/of the ad^ecant box and ^a ^oc^e\^ ^^l°cus of 
the shear center of the front spar and at the web of the rear s^ar. 

Each element is isolated and its ends are alternately assimed to 
be rigidly fixed. Calculations are then made of the forces and momen s 
aoMiS ^ the element at Its free and fixed ends neoeeaary to canae a 
nnlt^ertloal displacement vlthout rotation or twist ® ^ 

in the Plane of bending without vertical displacement or tUst, and a 
im'i t. twist without displacement or rotation in the plane ° "i^flu- 

These forces and moments at the ends of the elemen s ar 
ence coefficients/ 

displacements must eaual the external load on ^s_ 
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quantities in these equations in their final form is knovn as the "op- 
erations tahle." iiie solution of these Bimultaneous equations gives the 
displaceiaenbs at the ends of the assvmied elements for the given load 
condition. From these displacements and the influence coefficients the 
stresses at any point of the 8 t 3 ructure may he calculated. 

In the calc\ilation of the ini'lusnce coefficients the spaz’ elements 
vere first assumed to have shsai'— resistant vehs. Consequently ^ the mo- 
ment of inertia of the rear spar section was computed considering the 
four flange angles and the entire section of the wsh as fully effective 
in bending. In the front spar the semicircular leading-edge skin was 
also included in the effective cross-sectional area. Solution of the 
resulting equations gave emallei’ deflections for the front spar than 
those observed in the experiments, while in the case of the rear spar 
the agreement was satisfactory. For this reason the calculations were 
repeated on the assum.ption that the shear webs acted as diagonal— tension 
fields, and the angle of the diagonals was taken as 30 ° since the folds 
obseived in the test appeared to subtend approximately this angle. Cal- 
culations showed little effect of this angle on the displacements. This 
assumption resulted in good agreement between calculated and observed 
defiecfcions in the case of the front spar, but in slightly exaggerated 
calculated deflections in the case of the rear spar. It might be men- 
tioned here that the agreement between calculated and measured strains 
was generally better than that between calculated and observed deflec- 
tions . 

In the calculation of the influence coefficients of the interspar 
elements, namely, the boxes connecting the front spar with the rear ^ 
spar, the variation in the height of the cross section was duly consid- 
ered. 

CALCULATION OF THE INFLUENCE COEFFICIENTS 


Upon the beam of fig^ire 5 ^ shear force V and an end moment M 
are acting. It is kno^m tha-t the displacement v and the rotation m 
at the end can be calculated from the following formulas: 


= (l/3)VL°/EI + (1/2)ML^/EI 

(1) 

m = (l/2)VL^/EI + ML/EI 

(2) 


As was stated earlier, in the Southwell and in the PIBAL methods 
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onl;’’ one displaceniGnt Quantity at a tine 5 s considered different from 
zero. For iuetajice, it may te stipulated that the free end be displaced 
a unit distance in the vertical direction (downward) while the rotation 
is prevented by the constraints. Then with v = 1 eind n = 0 equa- 
tions (l) and (2) can be solved for M and V: 


M = -6EX/L2 

(3) 

V = 12EI/L® 

(4) 


These equations can be interpreted as expressing the values of the mo- 
ment M and the force V exerted by the constraints upon the free end 
of the bar when this end is d5splaced d.ownward and is prevented from 
rotating. By the definition given earlier, these aie influence coeffi— , 
cients . 

The influence coefficients correspcnding to a unit rotation and 
zero vertical displacement can be obtained in a cimilar manner, if in 
equations (l) and (2) v is set equal to zero and m equal to unity. 
The solution is 


M = 4EI/L , (5) 

V = -6EI/L2 (6) 

Influence coefficients are designated by two lower-case letters 
connected by an arc. The first letter refers to the force (or moment) 
exerted upon the bar, the ^cond to the displacement (or rotation) that 
caused it; consequently, mv is the moment caused by a unit vertical 
displacement, while ^ is the force caused by a unit vertical dis— 
p3.acement. In order to Indicate the points at which the force is act- 
ing and the displacement is taking place, two subscripts are used. The 
first subscript refers to the first lower-case letter in the symbol for 
the influence coefficient, and the second subscript to the second let- 
ter. With this convention equations (3) to (6) my be witten in the 
form 
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•> 

~ — <SEI/L^ 


Wgg - 12EI/L® 

> 

~ 4EI/L 




( 7 ) 


It may Be noted that mvgg = which equality is a con3eq\ience of 

Maxwell’s reciprocal theorem. 


The forces and moments exerted upon the fixed end of the Beam of 
figure 5 0 ’S rigid end reactions can Be calculated from the requirements 
of static equilibrium. They can Be expressed with the aid of the 
influence-coefficient notation, as follows: 


mv 


= mv-n-p = — 6EI/L*^ 




AB ~ -“^BB 




mm^ = mn^ = 2EI/L 




'^AB " ~^B 


(8) 


The 'Sign conventj on used in tliese equations considers as positive down- 
ward forces arw^ disp3.aceiiients as well as counterclockwise moments and 
rotations . 

As is wall known, equations (l) and (2) are vaJ.id only for Beams 


10 


NACA TN Ko. 1324 


of constant moment of inertia. 4>hen the moment of inertia varies, as in 
the interspar memhers, the deflections and rotations must be obtained by 
integration. Numerical examples are given in appendix A. The torsional 
rigidity of the rear spar is small because of the nature of the construc- 
tion. Advantage vas taken of this fact in the calculation of the 
influence coefficients by assuming this torsional rigidity eq.ual to zero. 
Accordingly, the rear spar cannot exert any end momenta upon the inter- 
spar members, so that it becomes unnecessary to consider the effect of a 
prescribed rotation at the rear spar in a vertical plane containing the 
interspar member. \Vhen the vertical displacement of the interspar mem- 
ber is prescribed at its end at the rear spar, the member deflects as a 
cantilever x^ith a concentrated load at its free end. blien a vertical 
displcicement or a rotation is stipulated at its front spar end, its rear 
spar end is considered as simply supported. 


In 8-ppendix B influence coefficients of Wagner beams are developed 
using strain-energy methods. With the notation of figure 6 the results 
are: 


“ y^BB ~ ^ (L^/Ah ) — (K/L) ] 


yy_ = 2E/(K - L^Ah^) 


'BB 


) (9a) 


ttg3 = (AEh72) [(1/L) - l 7(L® - h^AK)] 


vhere 


K = (8L/ht sin7d) + (4L73Ah^) + (L oot^'a/k) 
+ (2d^ tan^o/h^Ag) ') hi 


(9b) 


where A is the upper flange area, equal to that of the lower flange in 
this case. The influence coefficients at the fixed end can be calcu- 
lated from the following equations: 
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3^A3 = -2^BB V (10) 

^■^AB = 

The expreesioDS for yypB’ ^%’B-' y^BB given in equa- 

tions (9). The Influence coefficients of a twisted thin— walled bar can 
be calculated from Bredb*s formila 



( 11 ) 


where t is the relative angle of twist between the two ends of the bar 
and T is the torque. Tlie integration must be carried out arotind the 
entire perimeter. Setting t = 1 and solving for T yields 


T = 



4A? G 
f (ds/b) 


( 12 ) 


It is easy to see that the influence coefficients for the fixed and free 
ends differ only in sign. When the cross section varies^ the influence 
coefficients must be obtained by integration. For the interspar members 
this was done in appendix C. In view of the assumption that the torsion- 
al rigidity of the rear spar is zero, it is unnecessax’y to calculate 
influence coefficients corresponding to torsion of the rear spar. 

The calculation of the influence coefficients of the D-shape front 
spar under shear force and bending moment leads to some complications 
when the flat web develops diagonal tension. Obviously, part of the 
shear force is tramsmitted by the diagonal— tension field and the rest by 
the curved leading-edge covering. The leading-edge covering did not 
develop wrinkles in the experiments. The normal stress is carried by : 
the flanges and the curved sheet. The distribution of the shear force 
between the two elements was determined with the aid of Castigliano ’ s 
second theorem. Of the total shear force P the part teken by the flat 
web was P — P^. The total strain energ 7 stored in curved sheet. 
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flanges, weTj, and uprights, was calcu.lated and its differential coeffi- 
cient with respect to Px was set equal to zero. This yie3.ded an 
expression for Px in teiius of the load P, the end moment T, and 
the geometric and mechanical properties of the structure. 

Substitution of Px in the expression for the total strain energy 
and differentiation with respect to P and T led to the end deflec- 
tion and the end rotation, respectively, by vix’tue of Castigliano* s 
first theorem. The influence coefficients could then be calculated as 
was done in the case of the beam with sheai’— resistant web. Details of 
the calculations are given in appendix D. Nuaericai values of the in- 
fluence coefficients are collected in tables I to III. 


THE OPEEATIONS TASLE 


In the operations table are listed the forces and moments exerted 
upon each of the Joints of the structure as a result of linear and 
anguJ.ar displacements at the Joints. The symbols used to denote forces, 
moments as well as linear and angular displacements, and the sign con- 
vention adopted, are shown in figure 7. 

In order to explain how the entries in this table are calculated, 
let it be asstimsd that a pesitive unit rotation in the t— direction is 
undertaken at Joint 3- (See fig. 8.) For convenience the unit chosen 
is 10 ^ radian. The effect of this rotation upon all the Joints cf the 
structure must next be determined. However, as was stated earlier, 
whenever a displacement is undertaken, all the other possible displace- 
ments of the structure are assumed to be zero. In other words, the far 
ends of the members Joined at 3 are considered to be rigidly fixed: The 
effect of a rotation at Joint 3 is felt, therefore, only in members 
1-3, 3--4, and 3~5. 

It is apparent that the rotation stipulated will cause torsion in 
element 3~4. The moment that must be exerted upon this element to 
cause the prescribed rotation is the influence coefficient tt shown 
in figure 8. In addition, moments ttas are needed to rotate elements 
1—3 and 3~5 so that altogether a moment equal to 2 t ^33 + tt is re- 
quired. If the properties of bars 1—3 and 3— 5 were not the same, 
the moment tt 33 required to rotate each bar would be different. This 

is the case at Joint 5; since the length of element 3—5 is 24 inches 
while that of element 5—7 is 12 inches. As may be seen frcm figure 8 
the vertical forces that must act upon elements 1—3 and 3—5 at 
Joint 3 are equal in magnitude and opposite in sense. They add up to a 
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zero reoultaiit force. The forces and moments needed at the far ends of 
the e.lem0nts are also indicated in figure 8. The numerical values of 
all these iiiflaence coefficients are listed in tables I and III. 

Tlie operations table corresponding to the assumption of shear- 
resistant spar vebs is presented as table IV. Each of its columns 
contains the forces and moments at al3. the Joints of the structure 
corresponding to the displacement indicated at the top of the column. 

The location and the natui'e of the foi’ce (whether a force^ a moment in 
the n-direction;. or a moment in the t-direction) are shown at the 
left end of each row. In the operations table are listed the reactions 
of the forces and moments -wlrcli act upon the bars. Hence the entries 
in the operations table are the quantities calculated earlier multiplied 
by —1, In other words, the operations table contains the forces and 
moments exerted by the bars upon the assumed geometric constraints (the 
Joints) . 

In the co.luTnn headed t3 the entries in the first five rows are 

zeros. This corresponds to the fact that a rotation in the t— direction 
at Joint 3 has no effect upon Joint 2 because of the assumption of rigid 
end fixation at the far ends of the bars. (See figs. 4 and 8.) More- 
over, does not introduce a vertical force or a moment in the 

n-direction at Joint 4. The next item is T^, which is equal to 

112.86580 inch— pounds when the webs are assumed to be shes.r resistant. 
The positive sign corresponds to the sign convention of figure 7 *^0 

the fact that the operations table contains the forces and moments ex- 
erted upon the Joints (or constraints). Eight digits are given, since 
little additional work is involved in keeping a large number of digits 
when a calculating machine is used and since it it desirable to have a 
great accuracy for checking purposes. 

The next three rows refer to the effect of t3 upon joint 6. Three 
zeros are listed since t.3 does not affect joint 6. At joint 3^ T3 
is equal to -(2tt33 + tt) , as was explained earlier. The numerical 

value is ~407.425S0 inch— pounds, which can be checked easily with the 
aid of tables I and III. Again, as in this entire discussion of table 
IV, the numerical values are those corresponding to shear— resistant 
webs. In row 11 a zero is entered because the vertical forces at Joint 
3 cancel, as was sho^-m in the discussion of figure 8. In rows 12 to 15 
are entered the Influence coefficients for the fixed ends 1 and 5 mul- 
tiplied by —1. 

The entries in the other columns of the operations table are ob— 
tai.ned in a similar manner. In addition to the operations, table proper 
the so-called right— hand— side members are also listed in table IV . 

These represent the applied forces and moments acting on each of the 
Joints . 
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In the calculations, as well as in the experiments, four conditions 
of end fixation and loading -were conttidered. They are: 

Condition I.— Eoth spars are rigidlj- fixed to the test rig at their 

ends (joints 7 and 8 in fig. 4) and a 50C-pound load is acting downv/ard 
on the rear spar 5 inches outboard of Joint 1. llius, for this condition, 
the external ■vei'tical load on Joint 1 is 5'CO pounds and the ex';'.ernai 
moment in tlie t-di.rection is 2500 inch-pounds. The loads at all the 
other Joiiits ai’e zero. The operations tat.le corresponding to these end 
conditions is prese.rited in talale IV, and the externo-1 load.s are shown 
in the column, headed — R.H.S .j. 


Cond ition II .— The end fixation is the same as in condition I. A 

500--po\’nd load is acting dovmward on the rear spar at Joint 3. Hence 
the extemal load at 3 is 50C pounds, while the external moment at 3, 
as well as the external leads at a.ll the other Joints, is zero. The 
operations tab.le is contained in table IV and the external loads are 
listed in the column headed — E.H.S.jj, 

Condition III.— Tlie front spar is rigidly fixed at Joint 8 as be- 
fore, but the attachment of the rear spai- to the test r3g at Joint 7 is 
removed. The operations table is, of course, modified because of this 
change in the end conditions, but only’- the operations involving a ver- 
tical displacement y^ and a rotation tg at Joint 5 are affected. The 

quantities listed in table IV in the rows designated Tg and Yg as 
well as in the columns headed tg and yg must be replaced by the fol- 
lowing entries: 


Tg -260.145798 

Yg -9.20500 

The loading terms are the same as in condition I. 

Condition IV .— The end fixiation is the same as in condition III aud 
the loading is the same as in condition II. 


-9.20500 

-3.O81G933 
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Table P/ represents 2 sets of equations each containing 15 linear 
equations with 15 un.knowns corresponding to conditions I and II. For 
example, row 1 can be read as follows: 


-1325.907T32n2 - 54.5600T2ys + 39.926835^4 + 54.560OT2yi =0 (13) 


How 13 represents the following equation if condition I (corresponding 
to E.H.S.j) is considered: 

54.5600T2n2 + 2.31481^2 + 9.20500t3 + 0.767083373 

+ 9.20500ti - 3.081893571 = - 500 (14) 


The operations table based on the assumption of fully developed 
tension diagonal field action is presented as table V. It was con- 
structed in exactly the sarae manner as table IV, takirig, however, from 
tables I and III the influence coefficients corresponding to Wagner 
beams . 

The systems of linear equations were solved by ma.trix methods. Of these. 
Grout’s procedure appears to be most advantageous (reference 9). The 
displacement quantities obtained by solving the systems of linear equa- 
tions are listed in tables VI and VII for both assumptions of spar— web 
action. 


CALCULATION OF THE STRESSES 


ivhen the displacements and rotations are known at each joint, the 
forces and moments exerted upon the bars at the joints can be calculated 
easily with the a.id of the influence coefficients. I!hs shear force, the 
bending moment, and the torque at any point along the element can then 
be determined without difficulty from the laws of statics. 

As an example, the normal stress will be computed in the flange at 
the midpoint of element 3—5 for condition I on the assumption of 
shear— resistant webs. The displacement quantities needed are taken 
from table VI: 
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Jg = 0.39023 liiCh 

yg = 0.054731 inch 
t., = O.OI591S i-adlan 
tg = 0.0078959 radian 

The influence coefficients of element 3"“5 9-i’e rewritten from table 

^33 ~ 7670.8 pounds per inch 

y^.^g = -92,050 poimda per I’adian 
t^gg = 1,472,800 inch— pounds per radian 
yygg = --7670.8 pounds per inch 

^gg = -92,050 pounds per radian 
ttgs = 736,400 inch pounds per radian 
The shear force on element 3“5 nt Joint 3 iss 

P 3 = 73 (^ 33 ) + ■t3(yt33) + ys^yyas) + 45(yt35) 

The moment acting upon the element at Joint 3 is: 

I4g = ta^ttss) -f ysvtyss) + tsCttgg) + y5(ty3g) 


I: 


(15) 


(16) 


Substitution of the numerical values yields: 
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P3 = 381 po-unds 

M3 = —1624 inch-pounds 

Consoquently, the tending moment Mp at the midpoint of element 3—5 
is: 

:= 381 X 12 - 1624 2954 inch-pounds 


Since the moment of inertia of the rear spar is 0,84l6 inch"^ when the 
web is considered fully effective^ the maximum stress f^ in the flange 
is: 


fg = 2( 2954 )/( 0.4208) = 7020 psi 


DESCRIPTION OP THE EXEEEIMENTS 


In order to verify the theory, the aluminum -alloy model of a wing 
shown in figuire 1 was designed, constructed, and tested. The wing cont- 
prised two spars, each having a thin—sheet web and four angle— section 
flanges formed on a bending brake. Tb.3 curved leading— edge covering was 
attached to the flanges of the front spar so as to form a D-shape con- 
struction, The webs were stiffened by means of uprights riveted to them. 
Every third upright was part of a wing rib and was attached to the spar 
flanges as w^ell as to the web. Tension diagonals developed in the web 
under comparatively low loads, probably because the stiffening effect cf 
the shorter uprights was slight. 

The cun^ed portions of the flanges of the ribs were manufactui’ed on 
a hydraulic press and were obtained through the courtesy of Edo Aircraft 
Corporation. The rib webs as well as the upper and lower cover plates 
of the wing were also stiffened by light angle sections which, however, 
seemed to be sufficiently rigid to prevent buckling. At least no waves, 
diagonals, or wrinkles could be observed in the webs of the ribs and in 
the cover plates during the experiments. 

Three of the spaces between adjacent ribs were covered with skin, 
while the other three were left open to simulate cutouts. To the ends 
of the spars heavy machined steel fittings were riveted in order to pro- 
vide attachments to the test rig strong enough to represent rigid end 
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fixation. Of the two end attacliTiients the one at the front spar was al- 
ways bolted to the test rig, while the one at the rear spar was bolted in 
test conditions I and II, and free in conditions III and IV". 

Detai3.s of the construction and test arrangement may be seen in the 
photographs in figui’es’2 and 3. crder to provide for a suitable ap- 
plication of the concentrated loa.ds, stee3- channels wei*e riveted to both 
sides of the spar at locations corresponding to Joint 3 Q'-ad to a point 
5 inches outboard of Joint 1. (8ee fig. k.) Tlie external loads were 
applied to the steel channels by placing weights into a frame suspended 
from fittings attached to them. 

Strains were measured by means of Ba.ldwi.n Southwark SR— 4 type A— 1 
metaloctric strain gages cemented in pairs to the aft flanges (upper 
and lower) of the front spar and the forward flanges of the rear spar. 
Gages were located at the cutout sections of the wing at points A, B, 

C, D, E, and F in figure 4. The pairs of gages were connected in 
series in order to obtain the average value of the normal strain in each 
flange. An SR- 4 control box was used for measuring the strain, and a 
brass plug and tapered socket arrangement was employed for switching. 

The accuracy of the strain measurements was checked by tests made with 
a cantilever beam to which pairs of gages were cemented. The maximum 
error in strain was found to be about ±10 X 10 

The absolute values of the strains measirred in upper and lower 
flanges at the same location w^ere found to be very much the same. As an 
example, the strains observed in the rear spar in loading condition I 
are shown in figure 9« Figure 10 demonstrates that the variation with 
load of the average absolute value of the strain at any location was 
linear. Similar 1;^", the deflections increased proportionally to the 
loads, as may be seen in figure 11. 

The deflections were measu.red at points 1, 2, 3^ 4, 5^ Q^d o 
(fig. 4) by means of Ames dial gages placed on a sturd;^" steel frame 
rigidly attached to the test rig. 


COMPARISON OF THE RESULTS OF E3CPERIMENT AI® THEORY 


The final results of the experiments are presented in figures 12 to 
19, which contain the experimental curves of deflection and strain for 
the four conditions of loading. The values obtained by calculation^are 
also shown in the same figures reduced to corresjjond to a load of lyO 
pounds in order to facilitate comparison with the experimental results. 
From such a comparison the following conclusions can be drawn. 
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Calculated strains a^res veil with the experimental values for both 
spars in all the four conditions of loading when diagonal-tension action 
is asEuiaed in the spar webs. The oiily case in which there is appre- 
ciable disagreement is the strain in the rear spar for condition IV. 
However, in this case the strains are veiy small (at the point of great- 
est deviation 13 x 10"^ corresponding to about 13b-psi stress) and 
consequently herd to meas'ure and itninportaint for engineering purposes. 
Wlien the webs are ass\!med to be shear resistant, the agreement is .less 
satisfactory’’ for the front spar and not much different from that ob- 
tained by the diagonal— tension assumption in the case of the rear spar. 

As far as deflections are concerned, the diagonal— tension assump- 
tion gives good agreement for the front spar but not quite so good 
agreement for the rear spar for a.ll conditions of loading. The shear- 
i-e si stent- -web assumption is better for the rear spar and less satisfac- 
tory for the f?.'ont spar. 


C0NCI,US10KS 


The calculation of the strains in and the deflections of the spars 
of stressed— skin-type wings having cutouts can be caiTied out with a 
reasonable amount of work if the PIBAL method is used. Most of this 
■work may be done by little- ti’ained persorjnel. Time may be saved if use 
is made of the formu.las developed in this report for the influence co- 
efficients. Experiments with a model ^vlng gave satisfactory agreement 
with the theory for four different conditions of loading and end fixa- 
tion. 


Polytechnic Institute of Brooklyn, 

Brooklyn, E. Y., August 23, 1946. 
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APPEiroiX A 


DEFlECTIONS AWD EOTATIONS OF BEAM EIEMEKT 'v/ITH VARIABLE MOMENT OF INERTIA 


In calculating the vertical displacement at the aft end of the in— 
terspar element for a unit load at this end, the variable moment of 
inertia of the element (rib flanges, rib web, wing covering) between 
the rear spar and the front spar webs was calculated to be 

I^ = 0.0000l6x° + 0.00o24x^ + 0.2400x + 2 .U4l6 (Al) 


where x is measured from the rear spar. In the analysis based upon 
the asavanption of shear-resistant spar webs, the moment of inertia was 
assumed to vai*y hyperbolically between the front spar web and the shear 
center of the R-section, being Infinite at the latter point, and may be 
written as 


I^, = 3.57If/x' 


where x* is measured from the sheai-* center; 3.57 inches is the dis- 
tance between the shear center and the front spar web, and Ip is the 
moment of inertia of the interspar element at the front spar web . 

By snaking use of the unit— load method, the total displacement at 
the aft end for unit load at this end is expressed by 

20 3.57 

d = 1/E r dx/lx + (1/3.57 IfE) f X- (23.57 -x')^ dx* (A2) 

o o 

The first integration was performed numerically by Simpson’s rule, which 
is explained in reference 10. The numerical value of the vertical dis- 
placement is 


d = 43.20 X 10 ^ in. /unit load 
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To calculate the rotation at the forward end due to a imit moment 
in a vertical plena containing the interspar element, the unit— ].oad 
method 's again emploj'ed. For the as sumptions . of ahean-resistant spar 
webs and zero torsional rigidity of the irear spar, 


where cu is the angle of i-otation at the front spar in radians. The 
numerical value is 


It may be mentioned here that in the analysis for the assumption of 
Wa.;pie',’ beam action the interspar elements were assumed to have infinite 
rigidity between the front spar web and the shear center of the front 
spar. 

Becai’-se of the assumption of no torsional rigidity for the rear 
spar, for a typical interspar element 3—4 the influence coefficients 

S133, 1^33, 5^33, 1^34, 30x43, nri34, nn43 are zero. The remaining in- 
influence coefficients may' be obtajned from the values of d and co 
with the aid of formulas (A 4 ) . 


23.57 



3.57 

3.57 



(A 3 ) 


o 


CD = 0.07776 X 10 ~^ radian /unit moment 


^33 = -JJ 43 = 3^44 =■• -^34 = l/<i 


im44 = l/(X> 


(A 4 ) 


^44 = ny^44 = l/( 23 . 57 <u) 


yn34 = ny43 = - 23 . 57 /d 
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APPENDIX B 


VJAGNEPv BEAM INEniENCE COEFFICIENTS 


The total strain energy in the Vagner heam spar element shown in 
figure 6 may be expressed as 


U = (1/2E) jfjt htL + f (ft + ft^)^ At dx 

o 


f (f^ + <2^ + fy As 2 hi 

O 


(Bl) 


vhere t is the veh thickness. It may he stated that 

f 4 . = 2Y/(ht sin 2a) 

W \j 

ft = (Y/A^)[(x/h) - ( 1 / 2 ) cot a] 
f = -(Y,/A^)[(x/h) + ( 1 / 2 ) cot a] 
f^ = Yd tan. a/(hAg) 
ft =? T/(A^h) 

fc =-T/(Ach) 
o 




(B2) 


J 


The expressions for derived in the theory 

of tension field webs developed by E. Wagner and summarized by 
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Paul Kuim in reierenc© 11. In the exprsaslon for f . d is the distance 
■between npi-ights. 

On Euhstituting these expressions in equation (Bl) and integrating, 
the total strain energy hecomes 


1. ! , LY^ cot^ g 2YTL^ 

’ ht sin^2a 3Ah^ "■ aA ^2 


+ + Y"d^ tan^ a B hi 

Ah^ Agh^ 


(B3) 


In this equation A = Ag = A^. 

By virtue of Castigliano ’ s first theorem the partial derivative of 
U wtLth respect to Y yields the tota] vertical displacement d due 
to Y eind T, and the partial derivative of 'U vmth respect to T 
yields the total rotation t due to Y and T, d and E Being meas- 
ured at the free end of the element. 


“ = d = (I 'OE) r 8YL ^ 4L^Y LY cot^ g 2L^T 




$Y 


'ht(sin 2a) 3Ah‘' 


A Ah= 


2Yd^ tan^a Z hi f 
+ =• i 


> (B4) 


cU ^ t = YL^ + 2LT 
AEh^ 




Then with d = 1 and t = 0, as in the case of equations (l) and 
(2), equations (B4) are solved for T and Y. In the notation of the 
influence coefficients T = typg and Y = y^g. 
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In order to obtain the influence coefficients corresponding to a 
unit rotation and zero vertical displacement, equations (b 4) are solved 
d = 0 and t = 1. In the notation of the influence coefficients 

T = and Y = BbB = ^-yBB. The values of the influence coeffi- 

cients are given in equations (9) • 

APEEWDIX C 


IWFLUEWCS COEFFICISWT FOP TORSIOW OF BEAM FIEMEWT 
WITH VARIABLE MOMSHT OF IWERTIA 


Foi' the interspar members ^n this ana.lysis, the wall thickness of 
a cross section is constant, j> ds is replaced by P^. and A by A^^; 
P^ and A^ ai’e the perimeter and area of a cross section, respectively , 
at a distance x forward of the rear spar, for x less than 20 inches. 
In the notation of figure 20, 


Px = 4(0. lx + 8) 

Ax = 43(1 + 0.05x) 


(Cl) 


For X greater than 20 inches, again in the notation of figure 20, the 
p0xiTneter and area of a cross section of the member are given by 


P' = 2(2y + 12) 
A* = 24y 


( 02 ) 


Consequently, equation (12) becomes 

_ L o 


tt = 4Gb 


— dx + 4Gb 


o 


J 


r A’^ 


dx 


(03) 
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The TiTimerical valrie obtained le 


tt = 1,128658 X 10'^ in,— Ib/radian 


/iPPEM)IX D 


lEOET-SPAR IttmieHTCS COEFFICIEIJTS 


Figure 21 eho’/s a sketch of the cur\’’ed sheet and the flanges of the 
r— shape front spar, which carry the part Pi of the total vertical load 
P, and the end riomeut T. If f^ repi-esents the shear stress and f-j^ 

the bending stress, the strain energy' stored in the curved sheet is 


U 


curved sheet 


r> 

/ iidv 
■J 2G 
V 


r f= 


Y 


- dY 
2E 


(Bl) 


the integration being performed over the entire volume . Or, in the no— 
tatj on of f igiire 21, 


IT 


^curved sheet 


1 r f PiQ 


2G 


/ ( •— i 

-D \ lit / 


rtL dep 


L ^ 


1 f’ r f Kr cos cp \ 


+ “ 

2E 


o o 


V Ij 


/ tr dep dx 


where 
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<P 


f 


r^t cos 0 d0 = — + r^t sin cp 
2 


(D2) 


3, / , \2 

I, . 2?U. * 2A ( !i ) 

2 \ a / 

M = T + P3.X 


It will be noted that the flanges carry an additional normal stress 
due to ¥agner beam action. To calculate the strain energy stored in the 
Wagner beam a procedure similar to that given in appendix B is used. The 
strain energy in the web is 


U,., 


web 


U. 


flanges ^ 


4t 

2(P - Px) L 

Eht(sin 2a) 

in the f 

laoiges is 

L 

^ r / 

-,2 

f ( 

0 

, i'tlj 


L 

p 


r 


J ^ -C + ^’cx 


A dx 


(D3) 


(Di^) 


where f.^ is the normal stress in the tension flange due to V/agner beam 
action, f 4 . the normal stress in the tension flange due to bending re- 

suiting from (T + Pxx) of figure 21, and fc and f^ have corre- 
sponding meanings for the compression flange. From the discussion in 
appendix B it follows that 
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= 


I* ~ 1*1 f "x. cot a \ 
K - /' 




P — Pi / X cot a 
A 




-ti 


■Cl 


(T + Pix)h 
2I1 ‘ 


(T + Pix)h 

2 iT 


(P5) 




Substituting these values in equation (D 4 ) and integrating yields 


2. 3 


r2(P-Pi)"L" (P - Pi)^L cot^ a 

flanges ~ 1 p~p p 

2Ji >- 3A^h ' 


2(P -- Pl)L / qi 


AI 


r-i 


+ 

2 3 


+ — — f + TPiL + — — - ^ 
2ii^ ^ ^ 3 y 


(D6) 


The strain energy stored in the uprights of the Wagner beam is 


^uprights 


(P — Pi)^ d^ E hi tan^ a 
2Eh^Ag 


(D7) 


The total vertical displacement d and rotation t at the free 
end resulting from T and P are found from 
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au, . , 

^ _ tota ?, 
'~c>P 


"E" 


^TJ cTJ ^XJ ^XJ 

sheet ^ web ^ flanges ^ uprights ^ 

~ap ~ ap ' 8 p ap 

V 

/ 

^^^curved sheet ^^woh ^^^flanges ^ ^-'^uprights 


(D8) 


Ho^vever, it wi 3.1 he noted that the expression for contains Pi, 

the undotexminsd part of the total vertical load taken hy the curved 
sheet. Castigliano ’ 8 second theorem 


__ to tal 
"aPi 


(D 9 ) 


Men Pi is dete.ruiined from eqviation (D 9 ) the result substituted in 
the expression for Utotal> equations (d 8 ) yield d and t. Then pro- 
ceeding as in appendix B and substituting the numerical constants for 
the element tinder consideration yield 


~ 13,767.95 Ib/in, 


tyBB = -165,210.65 in.-lb/in. 


ttpB = 4,378,361.09 in. -lb /radian 


^AB = -(24ftBB + Ubb) 


= -413,305.439 in.-lb/radian 
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The nuTuerical constants used to obtain these results are: 

L = 24 in. 
h = 7.6 in 
L hi = 35 J in. 
t = 0.020 in. 
d = 4 in. 
a = 30 ° 

Aq = 0.06 sq in. 

A = 0.12 sq in. 
r = 4 in. 

= 5.4762 in.^ 

E = 10.5 X 10® psi 
G = 3.9 X 10 ® psi 
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TABLE I.-RSAE-6PAE lUFLIIEWCE COEFFICIENTS 


Infl\ience coefficients 

She ax'— re sis tant 
weh 

V/agner 

■beam 

Element l-S®’ 


7,670.833 

2,172.291 

tIfi, Yiti, T'S^i, Ylt3 

-92,050.000 

-26,067.488 

Titx, Tata 

1,472,800 

653,009.833 

Tits, Tsti 

736,400.000 

-27,390.169 

^a^s) ^3^1; Tijjj 

92,050.000 

26,067.488 

Y^-i, yI?3 

-7,670.833 

-2,172.291 

Element 5—7 

^sY 5 

61, 366,66 

6,135.8718 

TsJ'’5> Ygts, T7V5 

-368,200.000 

-36,815.231 

T?>5 

2,945,600 

901,291.383 

5 

—61, 360 . 6 b 

-6,135.8718 

T t 
7 5 

1,472,80c 

-459,508.615 

Yytg 

368,200 

-36,815.231 


^Fcr element 3-5 the Influence coefficients are 
identical with the corresponding ones given for element 1-3. 
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TABLE II FRONT-SPf\S INFINEKCE COEFFICIENTS 


1 

Influence coefficients 

Shear-resistant 

weh 

Wagner 

team 

Element 2—4 

YsYa, Y 4 Y 2 

62,066.211 

13,767.5543 

^ 2^4 

-744,794.531 

165,210.6517 


11,916,712,50 

4,378,361.087 

T'2t4, T 4 tg 

6 , 9 ?S, 3 ? 6.25 

-413,305.439 

T^4, E?' 4 , 4 j 

i 744,794.531 

165,210.651 

^ 472 .. ^274 

-62,066.211 

-13,767.5543 

Element 6-8 

_> 

Y V 
& e 

1 

496,529.69 

41,561.5189 

^'sJ '6 

-2,979,173.12 

-249,369.114 


23,833,425 

6,287,872.293 

"^ 6^6 

11,916,712.5 

-3,295,442.925 


-496,529.69 

-41,561.5189 

^8^6 

2 , 979 , 1 '^ 8.12 

249,369.114 


^or eleiaent 4-6 the influence coefficients are iden- 
tical with the corresponding ones given for element 2 - 4 . 
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TABIE III.- IIIFUJENCE COEFFICIEHTS FOB TEE INTERSPAR ELEMEOTB 


Influence coefficients 

She ax— re s i s t ant- 
wet assumption 

Wagner— beam 
assmaption 

Unit 

^rr = ^ff 

23,1^*10 

34,319.578 

Ib/in. 

ttff = ~btj.f = “ ~^fr 

1,128,658.0 

1,260,652.5 

in.— Ib/in. 

nn^^ = nn^P 

12,859,809 

17,613,632.5 

in.lb/radian 

= -y^f 

545,600.72 

732,132.04 

Ib/radian 


^Subscript f refers to front spar; subscript r refers to rear spar. 
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TABLE V.- OPERATIONS TABLE FOR THE ASSUMPTION OF WAGNER-BEAM SPAR ELD/ENTS. 
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TABLE VI 

SOIUTIOW OF MATRICES FOE TEE PSSUbmUON OF SHEAR-RESISTANT WEBS 


4 *1 

[Linear and Angu.Iar DispD.aoements times 10 J 


Condition 

I 

(500-lh load 
at end; Both 
spars fixed) 

II 

(500— lb load at 
joint 3; both 
spars fixed) 

III 

(500 "lb load 
at end; rear 
spar free) 

IV 

(500-lb 
load at 
Joint 3; 
rear 

spar free) 

ri2 radians 

lbO.24 

02.702 

310.93 

227.89 

,72 inches 

4619.7 

1720.1 

5,232.9 

2125.0 

t2 radians 

124.65 

39.849 

136.52 

47 . 401 

n4 radians 

84.315 

66.539 

233.83 

211.41 

y4 inches 

1919 - *+ 

792 . 4 o 

2,249,1 

1016.2 

t4 radians 

98.647 

37.092 

no . 36 

44,628 

ns radians 

14.342 

13.141 

147.57 

147 . 58 

ys inches 

239.49 

111.51 

307.13 

161.11 

t6 radians 

39.751 

18.435 

49.509 

25.174 

ts radians 

159.18 

66.529 

172.51 

74.191 

ys inches 

3902.3 

2388.0 

7,767.3 

6033.6 

ti radians 

179.62 

44.217 

191.69 

51.752 

yi inches 

8452.1 

3681.2 

12,619 

7508.5 

ts radians 

78.959 

54.528 

IOS.13 

68.645 

ys inches 

547.31 

401.39 

3,938.3 

3808.8 
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TABLE VII 

SOLUTION OF MATPICES FOR VfAGNEE-BEAM ASSUiVG'TION 


[Linear and Angt^lar Displacements times lo"^] 


Condition 

“r 

I 

(500-lh load 
at end; Doth 
spars fixed) 

II 

( 500-1 b load 
at Joint; loth 
spars fixed) 

III 

(500- It- load 
at end rear 
spar free) 

IV 

( 500-lb 
load at 
Joint; rear 
spar free) 

Us radians 

0.023300 

0 . 013069 

0.036220 

0.025316 

yg inches 

59731 

.22769 

.67828 

.28398 

t-2 radians 

.015446 

,0046967 

.016833 

.0059464 

radians 

.013163 

.012013 

.025987 

.021016 

inches 

.25218 

. 11187 

.29998 

.14094 

radians 

.012297 

.0046496 

.013664 

.0055411 

n^ radians 

.0033627 

.0032997 

.015063 

.015008 

Jq inches 

.032725 

. 01^46 

.048.969 

.030393 

tg radians 

.0050065 

.0023852 

.0061564 

.0030680 

t3 radians 

.00170180 

.0069400 

.018568 

.0078426 

ya inches 

. 56667 

, 40506 

.92556 

.64747 

radians 

.019828 

.0051996 

.021238 

.0060373 

Yx inches 

I 1.1635 

.54267 

1.5552 

.89538 

tg radians 

! .0081012 
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Figure 3.- Front view of wing model and test apparatus 
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Fig. 3 



Figure 3.- Rear view of wing model and test apparatus. 
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/.Z,3A.5,6, indicaie locatw/i of Ames ga<jes L =23.57 m first matrix 
Letters indicate location of strain gages 24.06 m second matrix 
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Figure 4.- Simplified structure assumed in analysis. 
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Figure 5.- Beam element for cailculation of influence coefficients. 



Figure 6.- Wagner beam element 
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FRONT SPAR 



REAR SPAR 


Figure 7.- Sign convention of forces, moments, linear and angular displacements. 

(All directions shown are positive; upper-case letters 
represent forces and moments; lower-case letters repre- 
sent vertical and angular displacements.) 
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Figure 8.- Effect of a unit rotation at joint 3 in the t-direction. 

(The influence coefficients shown represent the 
absolute values of the forces and moments exerted 
upon the structural elements by the constraints; 
the arrows indicate the directions in which they 
act . ) 
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riguTe 9,— Experimental etreilnB In upper and lower 
flanges, rear spar. Condition I. 
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Figure 10.- Variation of strain with load, from experiment, 
rear spar. Condition I. 

(Strains are averages of absolute values 
measured in upper and lower flanges.) 
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Figure 11.- Variation of deflections with load, from experiment, 
both spars. Condition I. 
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Figure 12.- Experimental and theoretical strains, front 
spar. Conditions I and II | load, 150 pounds. 
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Figure 13.- Experimental and theoretical strains, 
rear spar. Conditions I and II ; 
load, 150 pounds. 
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Figure 14.- Experimental and theoretical strains, 
front spar. Conditions III and IV; 
load, 150 pounds. 
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Figure 15.- Experimental and theoretical strains, 
rear spax. Conditions III eind IV; 
load, 150 pounds. 
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Figure 16.- Experimental and theoretical vertical 
displacements, front spar. Conditions 
I and II,* load, 150 pounds. 
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Figure 17.- Experimental and tbeoretical vertical 
dieplacements, rear spar. Conditions 
I and II; load, 150 pounds. 
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Figure 18.- Experimental and theoretical vertical 
dieplacements, front spar. Conditions 
III and IV J load, 150 pounds. 
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Figure 19.- Experimental and theoretical vertical 
displaLcementa , rear spar. Conditions 
III and IV; load, 150 pounds. 
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Figure 21.- Curved sheet and flanges of D-shape front spar. 

(Vertical loads and disolacements as well as moments and rotations 
are positive as shown.) 
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